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Abstract — We propose an expectation maximization (EM)- 
based algorithm for semi-blind channel estimation of reciprocal 
channels in amplify-and-forward (AF) two-way relay networks 
(TWRNs). By incorporating both data samples and pilots into 
the estimation, the proposed algorithm provides substantially 
higher accuracy than the conventional training-based approach. 
Furthermore, the proposed algorithm has a linear computational 
complexity per iteration and converges after a small number of 
iterations. 



I. Introduction 

Amplify-and-forward (AF) two-way relay networks 
(TWRNs) [1| require highly accurate information about 
the channel at the terminals, for both self-interference 
cancellation and coherent decoding. The problem of channel 
estimation for AF TWRNs has been considered in a number 
of works pl-p). Most works adopt a training -based approach 
that requires each terminal to transmit a pilot sequence known 
to the other terminal. Despite the robustness and simplicity 
of this approach, the transmission of known pilots consumes 
bandwidth resources, which undermines the spectral efficiency 
of TWRNs. 

One alternative to training-based estimation is semi-blind 
estimation |6| which, in addition to using pilots, also incorpo- 
rates the received data samples into the estimation. Semi-blind 
estimation requires a smaller number of pilots, which makes 
it more spectrally efficient than training-based estimation [6|. 

In this paper we derive a semi-blind expectation- 
maximization (EM)-based channel estimator for AF TRWNs, 
assuming reciprocal flat-fading channels. Each iteration of 
the proposed algorithm has a computational complexity that 
is linear in the number of data samples. Using simulations, 
we show that, even with a limited number of data samples, 
the EM algorithm provides substantially better accuracy than 
the training-based least-squares (LS) estimator. Only a small 
number of iterations is needed to achieve convergence. 

The rest of this work is organized as follows. In Section |H| 
we present the system model. In Section III we derive the 



EM algorithm. Simulation results are presented in Section IV 
Finally, our conclusions are discussed in Section [V] 

II. System Model 

We consider a half-duplex AF TWRN with two terminals 
71, 72 and a single relay 7Z operating in flat-fading channel 
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conditions. Each transmission period is divided into two 
phases: 

Phase 1: Each terminal transmits a block of L pilot sym- 
bols, followed by N data symbols. We denote by t\ 



[*i 



,ti L ] T and t 2 = [£21, •• ■ ,*2l] T the pilot symbol 



vectors of 71 and T2, and by S\ = [sn, . . . , sijv] and 
•S2 = [s2ij ■ • ■ , S2n] T the transmitted data symbol vectors 
of 71 and Ti, respectively. The data symbols S21, ■ ■ ■ ,S2N 
are equiprobably drawn from the set S = {£,1, ■ ■ ■ ,£,m} 
of size M. The average transmission powers of 71 and T2 
during pilot transmission are P\ and P2, respectively, i.e., 
t^ti = LP\ and if^ = LP%, For simplicity, we assume 
that each terminal uses the same power during pilot and data 
transmission, i.e., E{sf«i} = NP X and E {s§ s 2 } = NP 2 . 
The corresponding received pilot and data signal vectors at 7Z 
are rp = ht\ + gt2 + w, and rp> = hs\ + gs2 + n where 
u> and n are additive white Gaussian noise vectors with mean 
zero and covariance matriceqj a 2 Ip and a 2 In, respectively. 
Phase 2: The relay broadcasts the vectors Arp and Arp, 
in sequence, where A > is the amplification factor. The 
corresponding received pilot signal vector at 71 is 



y = Ah 2 tx + Ahgt 2 + Ahu + wi 
and the received data signal vector is 

z = Ah 2 S\ + Ahgs2 + Ahn + n\ 



(1) 



(2) 



where u>i and ri\ are also additive white Gaussian noise 
vectors with mean zero and covariance matrices a 2 Ip and 
<j 2 In, respectively. The average transmission power of the 
relay is maintained at P r over the long term by using the 
amplification factor A = \/ P + p r , tT a ■ We are interested in 

the estimation of the cascaded channel parameters a = h 2 
and b = hg, which are sufficient for detection. The channel 
coefficients remain fixed during the transmission of the L 
pilots and N data symbols. Moreover, the noise variance a 2 
is assumed to be known at 71. 



III. Proposed Channel Estimation Algorithm 

Let 6 = [a, b] T be the vector of unknown parameters that 
we wish to estimate. The observed vectors {y,z} represent 
the incomplete data and the data symbols S2 represent the 
hidden data. Hence, the complete data is {y,z,S2}, and the 

II denotes the L X L identity matrix. 



corresponding log-likelihood function (LLF) is we denote as Q(a;8^'), is given by 

£(y,z,s 2 ;9) = -N\ogM-(N + L)\og(na 2 (A 2 \a\ + 1))- Q( a; e W) = _iV logger 2 (A 2 |a| + 1)) ' 

CT(A|a| + 1)V ^ (\\Gy~AaGt 1 -AIt 2 +A 2 aXt 2 \\ 2 ^ 



G 2 a 2 {A 2 \a\ + l) 



N M \ 

An iteration of the EM algorithm, say the tth one, consists y y p (t) i ^ _ ^^ . _ ^ + j?^ 1 2 
of two steps. The first step, called the expectation step (E-step) ~"t ~J I 

consists of evaluating the expectation (10) 



where 

' N M \ N M 



Q(9;9^) = E{c(y,z,s 2 ;9)\y,z;9^} (4) 

of the LLF of the complete data, C(y, z, s 2 ; 9), with respect to \i=i ?=i / «=i j=i 

the conditional PMF /(s 2 |y, z; 9^') of the hidden data given 

the observations and the current estimate 9^' = [a^\ b^] T of jv M 

9. The second step of the EM algorithm is the maximization and X = V^ V"* p\ Ci s u + *2**i- 

step (M-step) which consists of maximizing the expectation i=i j=i 

Q(9:9^') with respect to 9 to obtain an updated estimate 

fl(*+i) ; P m order to maximize ( fl0[ ) with respect to a, we will maximize 

it first with respect to the phase <fi a = Za and then with respect 

9 {t+1) = argmaxQ(0-0 (t) ) (5) t0 l a l" Minimizin 8 (!H with res P ect t0 0a, we obtain 



9 

N M 



In our case, the E-step and M-step are as follows: </>i* +1) = ^ ~ Z ( E E $] ( GZl ~ v4X 0) * { A * X t,J ~ AGs i 
E-step: We hav^] ^i=ij=i 

-, +(Gy-^Tt 2 ) H U 2 A't 2 -ylG ! ti 

g(0;0«) = -(iV+i)log(7ra 2 (A 2 | a | + l)) 



Substituting cf> a ' into (jTUJ, we get the following function 



<T 2 (i4 2 |o| + 1) ,_, (11) 

N M 

y-Aati- Abt 2 \ \ 2 + ^ ^ /?g | z, - Aafe - A&& | 2 I that depends only on | a | : 

^ Q(| a |;0W) = -^log(^(A'|a| + l))- ^ l f- 2 ^ l ° l + ^ , 
where P\ ' is the posterior PMF of the ith data symbol during a [A \a\ + 1) 

the tth iteration, given by ^ ' 

where 

-| Zi -Aa<*' S H-A6«^| 2 TV M 



o- 2 (A2|o(*)|+l) ' " " SJI 1 iV JU 1 



•£ e ,2 (A 2 la (t) l+1) ^ "»'••« — -s«i - i=ii=i 

fc=l 

jv A/ 



M-step: We need to obtain the values a^ t+1 \ &(* +1 ) such that X^= ^^^^ ) |A 2 A , ^-AGs H | 2 + — 1| A 2 A't 2 -AGti| 
{ a (*+i), & (t+i)} = ar g/i max^g(0;0W). (8) 



0=[a : b] T " JV M 



Regarding b^ t+1 \ it can be easily verified that the value of b 

that maximizes Q(9;9^') for a given value of a is ^^ # ., 



(Gy - i4Jt a ) (^^^ - AGti) 



W M ,(*) 



E E Pi,j £j («i ^asii) + * 2 y ^a* 2 *2 Taking the derivative of ((T2jl w.r.t. |a|, we get 



M«) = ^^ ^ • (9) 



AEE /f j + At 2 ff t 2 dQ(|a|;gW) = iVA 2 ^H 2 + 2V>|-^ 2 ;7-2iy 

<=ii=i ' dial A 2 |a| + 1 cr 2 (A 2 |a| + l) 2 

(13) 
Substituting b (a) in place of b in <§), we obtain an updated Setting dQ(\aUe * ) = ^ wg ^ obtain ^ quadratic equation 
objective function that depends only on a. This function, which 

A 2 V>| 2 + (2V + NA 4 a 2 )\a\ + NA 2 a 2 - A 2 U - 2W = 0. 

2 We ignore the term TV log M since it has no impact on the solution. (14) 



Solving ( fl4| >, we finally get 

1 



1(4+1) - 



2A 2 V 



(2V + NA*a A )+ 



(2V + 7VA 4 ct 2 ) 2 - 4A 2 i>(7Vyl 2 CT 2 - A 2 U - 2W) 

'(15) 

As we can see from (|9j, (jTTJ and (JT3J, the computational 
complexity of each EM iteration is O(N), i.e., it is linear in 
the number of data samples. 

IV. Simulation Results 

In this section, we investigate through simulations the MSE 
performance of the derived EM algorithm. We model h and 
g as independent complex Gaussian RVs with mean zero and 
variance 1. We assume that Pi = P-2 = P r and average our 
results over 100 independent channel realizations. The data 
symbols are generated from square QAM constellations. We 
consider the modulation orders M — 4, 16, 64. The pilot 
vectors t\ and £2 are obtained using M = 4 and chosen to be 
orthogonal. In our plots, we consider the total MSE, which is 
the sum of the MSE for the estimation of a and b. 

In Fig. [T] we plot the MSE performance of the derived 
semi-blind EM algorithm versus SNR for L = 8 and N = 32, 
where the SNR is defined as 10 log 5|. The channel estimates 
are obtained after 4 iterations of the EM algorithm. For 
comparison, we also plot the MSE of the LS estimator that 
only uses the pilot samples. As we can see from Fig. [T] the 
EM algorithm provides substantially higher accuracy than the 
LS estimator. Moreover, the gain in accuracy depends on the 
modulation order: the lower the modulation order the higher 
the gain. 

We finally consider the convergence of the EM algorithm. 
In Fig. [2] we plot the MSE of the EM algorithm versus the 
number of iterations for N = 32 and N = 100, assuming 8 
pilots and an SNR of 15dB. Fig. [2] shows that the number of 
iterations needed for convergence is small for all modulation 
orders. In all cases, convergence is achieved within at most 12 
iterations (as few as 4 iterations are sufficient in some cases). 
Convergence becomes slightly slower as the modulation order 
increases and as the number of data samples increases. 

V. Conclusions 

In this work, we derived the EM algorithm for semi- 
blind channel estimation in AF TWRNs assuming reciprocal 
flat-fading channel conditions. The resulting computational 
complexity of the EM steps is linear in the number of data 
samples. Using simulations, we showed that, even with a 
limited number of data symbols, the derived EM algorithm 
provides a substantial improvement in accuracy over the pilot- 
based LS estimator. Moreover, it requires only a small number 
of iterations to converge. 
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Fig. 1 . MSE performance of the EM algorithm for reciprocal channels along 
with the corresponding semi-blind CRBs plotted versus SNR (JV = 32 and 
L = 8, 4 EM iterations). 
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Fig. 2. MSE performance of the EM algorithm for reciprocal channels plotted 
versus the number of EM iterations (TV = 32, 100, L = 8, SNR 15dB). 
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